Conservative Estimates of the Mass of the Neutrino from Cosmology 
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A range of experimental results point to the existence of a massive neutrino. The recent high 
precision measurements of the cosmic microwave background and the large scale surveys of galaxies 
can be used to place an upper bound on this mass. In this paper we perform a thorough analysis of 
all assumptions that go into obtaining a credible limit on ^2 m v In particular we explore the impact 
of extending parameter space beyond the current standard cosmological model, the importance of 
priors and the uncertainties due to biasing in large scale structure. We find that the mass constraints 
are independent of the choice of parameterization as well as the inclusion of spatial curvature. The 
results of including the possibility of dark energy and tensors perturbations are shown to depend 
critically on the data sets used. The difference between an upper bound of 2.2 eV, assuming generic 
initial conditions, and an upper bound of 0.63 eV, assuming adiabaticity and a galaxy bias of 1, 
demonstrate the dependence of such a constraint on the assumptions in the analysis. 



I. INTRODUCTION 

The neutrino is an integral component of the stan- 
dard model of particle physics. Until recently it was 
assumed to be massless. With the new advances in 
non-accelerator particle physics there is now definitive 
evidence that this cannot be so: neutrinos must have 
mass. The first signatures for their masses were observed 
flavour oscillations in atmospheric and solar neutrinos 
and have later been verified in accelerator and nuclear 
reactor sources [l[ ■ When these observations are inter- 
preted within the 3-neutrino type scenario of the stan- 
dard model, their masses are required to be of the order 
of the measured mass differences, leading to two possible 
mass scales; 5m a t m — 3 x 10~ 3 eV 2 and , in conjunction 
with the nuclear reactor experiment results from Kam- 
LAND, Sm so i ~ 5 x 1CP 5 eV 2 . These mass scales can be 
accommodated by a model in which at least two eigen- 
states have mass. Alternatively if their true mass exceed 
O.leV then these results indicate that all three species 
are nearly degenerate. Tritium decay measurements have 
also been able to place an upper limit on the electron neu- 
trino mass of 2.3 eV at the 95% confidence level (CL) 
Although difficult to reconcile with the rest of the data, it 
is interesting to note the result from the Los Alamos Liq- 
uid Scintillator Neutrino Detector (LSND) which implies 
an lower limit of m v > 0.4 eV Q. Given our certainty 
that they exist, these results make the neutrino one of 
the most convincing dark matter candidates 

Although it is unlikely that these massive neutrinos are 
the dark matter it appears conceivable that they may still 
affect the growth of density perturbations in a measur- 
able way. Neutrinos with a mass less than 2 eV are still 
relativistic when entering the horizon for scales of k = 0.1 
h Mpc -1 and are quasi-relativistic at recombination. 
Therefore they cannot be treated as a non-relativistic 
component of the CMB and are not entirely degenerate 
with the other relativistic components In a matter 
dominated universe, a modicum of massive neutrinos will 
free-stream on scales of clusters and galaxies and there- 
fore suppress the rate of growth of density perturbations 



from being proportional to the scale factor a to being 
proportional to a 1 '" where e = |[1 - (1 - 240^/25) 1/2 ] , 
where the neutrino mass, ^2 m v is related through the 
expression: il„ = m„/(93.15/i 2 ). Here the Hubble 
constant today is Hq = 100/ikm s _1 Mpc -1 . This rela- 
tion can be applied even in the event that the neutrino 
species are found to be non-degenerate [3| ■ This effect on 
the growth of structure supplies us with a useful method 
for constraining f2„ and as a result, ^ m v . By measuring 
the amplitude of clustering on large scales (above the free 
streaming scale) and comparing it to the level of cluster- 
ing on small scales (below the free streaming scale) it 
is possible to tease out the level of damping due to the 
neutrinos. The amount of clustering on large scales is 
well constrained by measurements of fluctuations in the 
Cosmic Microwave Background (CMB) which map out 
the density perturbations on scales up to the horizon. 
Surveys of galaxies allow us to pin down the amount of 
clustering on small scales. Combining the two allows us 
to place a constraint on il„. 

This approach has been applied extensively over the 
last decade. In 0], using a combination of data 
from the Wilkinson Microwave Anisotropy Probe 3-year 
(WMAP3) 0, 2dF Galaxy Redshift Survey [|, Sloan 
Digital Sky Survey (SDSS) @ and the supernova (SN) 
[Til [TlT |. the WMAP team placed an upper limit of 
^m„ < 0.66 eV on the neutrino mass when the SDSS 
bias constraint f{| is included. Seljak [l2j recently re- 
ported an huge improvement on this when the same 
dataset is supplemented with the baryon acoustic oscil- 
lation (BAO) [HI and Lyman-a [l4[ constraints. The 
Lyman-cv forest provides information about the matter 
power spectrum on small scales (where neutrinos sup- 
press power) and at high redshift (z = 2 — 4) where 
the non-linear evolution is less significant. The Lyman- 
a data prefers a higher normalization and is the primary 
source of the strong upper bound on the neutrino mass of 
J2 m v < 0.17 eV at 95% CL in the absence of the SDSS 
bias constraint. This constraint is found using the recent 
measurement of the position of a peak in the galaxy corre- 
lation function at ~ 140 Mpc h -1 arising from the baryon 
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acoustic oscillations once they have decoupled from pho- 
tons. Given its dependence on Q, m and the hubble pa- 
rameter H (z) and the degeneracy of neutrino mass with 
both these parameters the use of the SN data to fix Q. m 
in conjunction with the information about £Ia provided 
from the BAO result proves to be powerful. 

The use of cosmological observations can, in princi- 
ple, supply us with limits on the neutrino mass which 
are comparative with experimental bounds. Yet there 
is a valid concern that cosmological constraints depend 
largely on the chosen cosmological model and accompa- 
nying assumptions/priors . The impact of certain as- 
sumptions that go into such an analysis has been consid- 
ered by other authors [HI, [l6[ , however a detailed review 
of all relevant factors that may be getting in the way 
of a truly robust constraint. In this paper we wish to 
elucidate how robust the neutrino constraints are to the 
assumed model. The structure of this letter is as fol- 
lows. We first establish the details of our the analysis 
and determine the way in which choices of prior prob- 
ability distributions for the neutrino content affect our 
results within the context of a ACDM universe. We then 
extend the model to include the possibilites of dark en- 
ergy, curvature, tensors, spectral running and more gen- 
eral initial conditions. We then consider how knowledge 
of cosmological parameters such as the galaxy bias can 
be harnessed to improve the neutrino mass limits. 



II. THE APPROACH 

Neutrinos will impact cosmological observables in a va- 
riety of ways. Relativistic neutrinos have a marginal 
effect by increasing the radiation density before decou- 
pling thereby impinging on the acoustic peak morphol- 
ogy. Massive neutrinos that become non-relativistic prior 
to recombination suppress the early integrated Sachs- 
wolf effect Their signature is a modification of the height 
and position of the first peak. Furthermore, the free- 
streaming of massive neutrinos induces a more rapid de- 
cay of the gravitational potential, fueling the acoustic 
oscillations within this free-streaming scale. This mod- 
ifies the heights of the second and third peaks relative 
to the first which can not be reconciled in the context of 
the shift in peak position and facilitates the constraint. 
As shown in 17[ it is thus possible to derive limits on 
neutrino mass from the CMB alone, however these above 
mentioned degeneracies weaken the power of the data and 
the constraints come out to be Yl m u < 2.0 eV for flat 
power-law ACDM. The merits of this are that the limit 
suffers less from systematic errors and is robust since 
it is derived from the single experiment with one set of 
systematics. The most effective use of the CMB data in 
pining down ^2 m u is however in the normalization of the 
large scale galaxy clustering power. The primary effect of 
the neutrino mass on the matter power spectrum P(k) is 
to reduce power on scales smaller than the free-streaming 
scale. Using both LSS and CMB data concurrently has 



been identified as the means of arriving at reliable and 
competative constraints. 

The empirical relationship between the morphology of 
the supernova type la light curves and their intrinsic lu- 
minosities facilitates their use as standard candles. These 
datasets [13, HU have become key in determining the ex- 
pansion rate of the universe and thus in testing models 
of dark energy. However the physical mechanism that 
forms the basis for this relation has yet to be established 
and given the enormous impact of these measurements 
we choose not to include them in our standard data set 
but will incorporate the SN results taken from [ll[ at 
points for the sake of comparison. 

We approach the Lyman-a measurements |14| with 
even more caution given the as yet still unresolved sys- 
tematics that still plague these data sets. When re- 
analysed b y 11211 and [18| discrepancies arose, with the 
authors of [1 8f | finding a lower normalization than the 
primary result in [lij . Since the Lyman-a data has been 
shown to tighten the bound on neutrino mass consider- 
ably we choose not to include it in our analysis. The 
impact of the inclusion of these different datasets on the 
neutrino mass limits has been comprehensively reviewed 
in [19( | and found to be significant. 

The practicalities of the method are as follows. CMB 
polarization and temperature power spectra as well as 
the matter power spectra are computed using the CAMB 
package [20]. We compare the spectra computed from 
the sample model to our fiducial dataset, comprising 
of the CMB temperature anisotropy measurements from 
WMAP 3-year [Tj and a combination of small scale (high 
i) CMB data from ACBAR, BOOMERANG, CBI and 
VSA [2l|, complemented with the galaxy power spec- 
trum derived from the SDSS [§]. We additionally use the 
constraint on the baryon density today of fl^h 2 from Big 
Bang neucleosynthesis (BBN) [H, [23|. In [I| tension 
between the WMAP 3-year data and the BBN measure- 
ment was identified, however this prior is necessary to 
exclude wild high fl^h 2 cosmological models which are 
favoured when isocurvature is allowed. 

A likelihood analysis using the likelihood function in 
[25l ] is then performed in order to compare the spec- 
tra generated from the models with the data. We use 
a Monte Carlo Markov chain method that invokes a 
Metropolis algorithm as described in [2(| H3] to explore 
the resulting likelihood distribution in parameter space 
efficiently. 

We take the concordance model as our starting point: 
a spatially flat Universe with nearly scale invariant adi- 
abatic fluctuations dominated by cosmological constant. 
A Fisher matrix analysis fl5j reveals strong degeneracies 
between the neutrino density f2„, the galaxy bias b and 
the matter density f2 m = + O^- These parameters im- 
pact matter power spectrum strongly on scales smaller 
than the free-streaming scale, while their effects on the 
CMB are similar in their subtlety (in the case of the bias, 
none at all). Neutrinos with masses of order eV become 
non-relativitistic only after decoupling, when the evolu- 
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tion determining the shape of the CMB spectrum has 
already taken place [J]. 

The bias parameter shares an effect on the normaliza- 
tion of the power spectrum on all scales so an increase 
in m f could be partly compensated for by a decrease 
in b. Increasing the matter content Q m brings matter 
domination forward in time, resulting in less suppression 
on small scales (i.e. shallower gradient at larger k [28|) 
and a smaller horizon at matter-radiation equality [291 ] - 
It then makes sense that the power suppression at larger 
k caused by the massive neutrinos must be corrected for 
this opposing effect; ~ —8^. A change in Q m is 
equivalent to a change in fl\ and the degeneracy between 
this parameter and the neutrino mass. 

Adding in 3 types of massive neutrinos degenerate in 
mass (we in fact use N = 3.04 as predicted by the- 
ory Q ) , the cosmological parameters that sufficiently de- 
scribe this scenario are the physical baryon density cjj,, 
the physical cold dark matter (CDM) density u)d, the 
neutrino density relative to CDM /„ = uj v /ujd, the frac- 
tional density of cosmological constant Q\, the galaxy 
bias b (assumed to be a constant based on theoretical 
reasoning in [3(|), the amplitude and scalar spectral in- 
dex of the primordial fluctuation spectrum A s and n s 
respectively. 

The parameterization used by the WMAP team f v is 
the ratio of the density of this hot dark matter compo- 
nent to the CDM. Given the degeneracy between uj v and 
Qd the ratio of the two is likely to be a sensitive param- 
eter. There are however alternatives; in other work the 
parameterization is chosen to be F u , the fraction of the 
total matter density fl m = flb + ^d + ^v that the neutri- 
nos comprise. In other cases one simply works with the 
total mass, ^m„, of all 3 species. It seems worthwhile 
to check the effect of the neutrino parameterization by 
repeating the same analysis using Y m v and F v = 
as the variables. If seemingly benign flat priors are im- 
posed different parameterizations the effective priors on 
the quantity being constrained, in this case ^m„, will 
diverge. Taking the priors shown in table [ill which are 
typically placed on the various neutrino variables in such 
work, we find that the upper limits on Y fn u at the 95% 
CL (using the CMB and LSS data alone) are in agree- 
ment for all 3 parameterizations. In the region of Y, Tn u 
space in which the data falls (< 4 eV), these priors are all 
effectively equivalent top-hat and thus have no impact. 
If we add the SN data [11| which amounts to imposing 
the constraint Q, m — 0.263 ±0.042, the effective posterior 
probabilities look quite different in this region of Y, m v 
space, but again do not lead to major discrepencies in the 
mass constraints. The situation does not change when 
a(Q m ) is reduced to 0.01 as predicted for future experi- 
ments such as SNAP. This is reasonable given that FIm 
and fib are well constrained by the current data. We con- 
clude the use of different subsets of the current data does 
not affect the independence of parameterization on the 
neutrino mass limits and thus continue using /„ through- 
out the rest of the analysis. 



Parameterization Prior m " f° r N v = 3.04 

95% CL 

f v = fi„/fid - 1 1.26 eV 

F v = Q V /Q M - 1 1.27 eV 

J2m„ - 5 eV 1.29 eV 



TABLE I: Constraints on neutrino mass when different pa- 
rameterizations are used. 

We now explore deviations from the standard model, 
specifically spatial curvature characterized by £Ik, the 
presence of dark energy, the presence of tensors, the run- 
ning of the spectral index and isocurvature contributions 
to the initial fluctuation spectrum. 



III. EXTENDING PARAMETER SPACE 
BEYOND ACDM 

Despite the vast amount of data now in support of 
the simple concordance picture, more complicated sce- 
narios are possible in certain paradigms. In order to be 
certain that the least restrictive mass limit is obtained, 
the model-dependence of these constraints should be ex- 
plored. 



A. Dark energy? 

In recent work [f| [U [HJ a degeneracy between the 
equation of state of dark energy w and the neutrino mass 
was identified. Massive neutrinos contribute to the en- 
ergy density of the matter component in the universe and 
thus can potentially alter the epoch of matter-radiation 
equality as well as the diameter angular distance to the 
last scattering surface. We wish to clarify this result. If 
w is allowed to vary, the effect of larger neutrino masses 
can be mimicked in the power spectrum by more negative 
values of w instead of larger £l m values which become in- 
compatible with the SN data [3l[ . If we make no assump- 
tion about the cosmological constant and include a con- 
stant w in the parameter space, the resulting confidence 
contours for the CMB and LSS dataset (disregard the 
SN results) show no signs of a relationship between these 
parameters. The inclusion of the SN data [ll| tightens 
the correlation between w and fl\ and hence links more 
strongly to Y m v, giving rise to the degeneracy between 
that has been pointed out. This contrast is demonstrated 
in figure ([1]) . The investigation of the correlation between 
neutrino mass and dark energy has been extended in [33j 
to include dynamical models in which the equation of 
state is characterized by a smooth time- varying function. 
Dynamical dark energy will impact the time-evolution of 
the gravitational potential wells and thus contribute to 
late time integrated Sachs Wolf (ISW) effect which will 
in turn manifest in the CMB power spectrum. A strong 
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FIG. 1: Extending the model beyond ACDM. From top right, 
bottom left to right: 95 and 68 % contours in the (^m^, CIk), 
(X m "> Qs ) an d (X^ m f r ) planes. Top left: The confidence 
intervals in the (Y^ra^, w) planes for our CMB+LSS dataset 
(blue) and for CMB+LSS+SN dataset (red). 



anti-correlation between ^ m„ and the time-dependent 
part of w(z) was found which accounts for the substantial 
weakening of the upper bound to vn v < 2.8 (depend- 
ing on inclusion/exclusion of dark energy perturbations) 
for the CMB+LSS+SN dataset. As our understanding of 
dark energy improves, we may have more reason to use 
more complex models of this energy density component 
into neutrino mass analyses. 



B. Spatial Curvature? 

The analysis so far has been performed under the as- 
sumption of fiat space. We now consider the possibility 
of non-flat geometry and include the curvature density 
VLk as a parameter. When we admit spatial curvature, a 
slightly closed universe is favourable with a best-fit model 
of Qk = —0.05. A tighter constraint on the neutrino 
mass of 1.17cV is also obtained. The CMB feels the im- 
pact as the inclusion of VLk which dictates the position 
of the first acoustic peak. Since Vl m h? also affects all 
the peak positions, the favouring of an open Universe by 
the data leads to a higher distribution in Q, m h 2 . But 
because of the increased freedom on the right-hand side 
of Vt m + fi r = 1 — D,k which holds for early times, a 
higher fl rn is countered by a more negetive £Ik , produc- 
ing a more gentle rise in the matter-radiation ratio and 
thus disfavouring the high ^ m„ models. The results are 
however consistent with flat space and massless neutri- 
nos within 2-ct. The confidence contour in figure |T]) rules 
out any degeneracy between the neutrino mass and Qk- 
This parameter does not have notable effects on the al- 
lowable models and can be confidently ignored from the 



parameter space. 



C. Running spectral index? 

Inflationary cosmological scenarios seem to provide the 
most convincing explanations for the initial conditions 
for structure formation [13. Although a scale-invariant 
spectrum of primordial fluctuations is regarded as an in- 
trinsic accompaniment of these models, this picture is 
only valid within certain cosmologies, such as pure de 
Sitte [34j . There have been claims of strong observational 
indications from the CMB, LSS and Ly-q data that the 
spectrum is a function of k (e.g. see [35j), however the 
fit of combinations of other datasets have been shown 
to deteriorate when running is allowed Q. Introducing 
running also means that the constraints on other param- 
eters are weakened and that the contribution of gravity 
waves can no longer be probed using the CMB data at 
all. In particular, there exists a strong degeneracy be- 
tween Q m and variations in the spectral slope on the 
scales probed by the CMB (l7l ? ].[36j. A spectral index 
which decreases on small scales (negative running) would 
limit the amount of small-scale damping that could be 
tolerated from massive neutrinos. Given that significant 
running of the spectral index is still under debate and 
has implications on neutrino mass constraints, it seems 
worthwhile to explore. In recent work [38| the neutrino 
mass constraints were found to be lowered in the pres- 
ence of a scale-variant n when the CMB, LSS and SN 
data was used. The SN data prefers a lower value of fl m 
which means that the effect of increased m„ can not 
be absorbed by an slight increase in the matter density. 
We consider a running spectral index in the presence of 
massive neutrinos without the admission of tensors using 
our fiducial dataset which does not include the SN re- 
sults. In this analysis we assume the form of the power 
spectrum to be 



d lnA|(fc) 
d lnfc 



n(fc ) 



d n s 
d ln(Jfe) 



ln(fc/fc ). (1) 



and parameterize the running of the spectral using a s — 
dn s /d\nk. The pivot scale is k = 0.05 Mpc -1 . We 
find the constraint on J2 m f opens up to 1.66 eV (95% 
C.L.) while the amount of running allowed comes out 
as a s = —0.062 ± 0.02 in agreement with the findings 
of the WMAP team using the 3-year data. This result, 
although at odds with the findings in [38[, is expected 
given the signficantly higher values of Q m that are pre- 
ferred in the presence of a s than in the case of constant 
n s . Without the SN data to restrict the matter den- 
sity, small-scale structure suppression can be offset by 
an increase in the matter density. Figure JTJ) shows the 
absence of a degeneracy between a s and ^2 m v- 
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D. Primordial tensor perturbations? 

Temperature anisotropics are sourced not only by den- 
sity fluctuations (scalar modes) but by tensor modes on 
small angular scales. Inflationary models predict their 
production by gravitational waves and evolution inde- 
pendently of scalar modes, leading to an uncorrelated 
power spectrum [3§|. On scales of the Hubble radius 
such modes interfere with the photon propagation along 
the line of sight and in so doing, induce extra anisotropy 
predominately on large angular scales. In the CMB, the 
presence of tensor perturbations is felt on small angular 
scales with increased low-£ power as well as a lower am- 
plitude of density fluctuations A s . Since A s affects the 
height of the first acoustic peak we can use the LSS data 
to normalize the CMB power spectrum thereby finding 
the best fit to the peak and the large scale part and in so 
doing can constrain the amplitude of the tensor contribu- 
tion, At , and A s simultaneously. Different models predict 
varying amounts of tensor perturbations and hence it is 
parameterized in terms of its ratio to the scalar mode 
anisotropies. Given the degeneracy between A t of these 
metric perturbations with w m , Wf, and n s [? ] it is impor- 
tant to check whether a higher contribution r is tolerated 
in the presence of massive neutrinos. We use the same 
convention as given in [40l | where r is the ratio of the pri- 
mordial amplitudes of the tensor and scalar fluctuations 

_ Ptensjk*) (l) s 

p , (it ) { ' 

1 scalar V™* J 

Here k* = 0.002 Mpc . This parameter relates to the 
spectral index of the tensor modes in the following way; 

n t = -r/8. (3) 

We do not allow for the running of the spectral in- 
dices which would reduce the amount of tensors that 
can fit the data. The upper limit on the neutrino mass 
degrades slightly to ^ TO ^ = 1-38 eV with the intro- 
duction of tensors with their fraction characterized by 
r = O.IOIq oi' Bearing in mind that the current theoreti- 
cal models which support the presence of tensors are few, 
it is reasonable to conclude from this result that exclud- 
ing tensor perturbations from the cosmological model will 
not effect the robustness of neutrino mass constraints. 



E. Mixed Initial Conditions? 

Up to this point we have assumed that the initial con- 
ditions in the universe are purely adiabatic. When this 
assumption is relaxed, 4 isocurvature regular modes are 
assumed to arise, namely; CDM (CI), baryon (BI), neu- 
trino density (NID) and neutrino velocity (NIV) (where 
the momentum densities of the neutrino and photon- 
baryon fluids are assumed to cancel exactly |4JJ). It 
was proposed in [42j that the CMB and baryon isocurva- 
ture modes can not be differentiated observationally and 



confirmed by the identical (up to an unimportant multi- 
plicative constant) correlation matrix elements [43j. We 
consider the simultaneous mixing of all distinct modes so 
there will be 10 components of the initial power spectrum 
matrix; 

< Ai(k)A*(k') >= P t3 (k)S 3 (k - k'), (4) 

where — 1,2,3,4,5 ) label the modes and their am- 
plitudes Ai. Assuming that the initial power spectrum 
of each mode is a smoothly varying function of k and 
given by Pjj (ft) = Aijk 7ls , where each correlation is char- 
acterized by an amplitude Aij and a single scalar spec- 
tral index (the auto-correlated adiabatic mode is already 
characterized by the parameter A s ). Following [45j it 
holds that 

Aij oc Zij (5) 

where 

N 

trace (zz T ) = z\j = 1 (6) 

i J = l 

So to characterize the all distinct isocurvature modes 
we introduce 9 additional parameters, namely the co- 
efficients z^ which give the relative contributions of each 
spectra to the normalization. This parameterization is 
advantageous because in our case where n aa < — rii SO , it is 
independent of the value of the pivot scale fco- 

1. Adiabatic + CDM isocurvature 

We initially consider the simplest extension of the adi- 
abatic regime, which is a sing le isocurvature mode mixed 
with the adiabatic mode [44[ . Including only the CI and 
AD modes we arrive at the same constraints as in the 
adiabatic case of 1.27 eV. We do not find any inter- 
dependence between the parameters characterizing the 
CI modes and the neutrino mass. 

2. AD + CI + NID 

The neutrino isocurvature modes are generally diffi- 
cult to motivate theoretically. Vanishing neutrino chem- 
ical potentials mean that there are the same numbers of 
neutrinos as anti-neutrinos produced in each generation. 
The implication is that spatial differences in the densities 
of neutrinos and photons will be naturally erased via the 
processes involving v v annihiliation. The netrino modes 
can thus only be generated at temperatures as low as a 
few MeV where these processes have ceased. Avoiding 
these issues requires the introduction of non-zero chemi- 
cal potentials. We find that when the NID isocurvature 
mode is admitted, significantly higher neutrino masses 
with an upper limit 2.4 eV at the 95% CL are permis- 
sible when compared to our adiabatic constraint of 1.3 
eV. 
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FIG. 2: ID marginalized parameter distributions for adiabatic 
only (solid line) and a correlated mixture of adiabatic and 
isocurvature initial conditions (dotted line). 



Additional parameters 



X)m„ for N v = 3.04 
95% CL 



Spatial Curvature Q.k 


1.17 


eV 


Dark energy w (w — constant) 


1.18 


eV 


Tensors At 


1.38 


eV 


Running spectral index n run 


1.66 


eV 


Isocurvature (all modes) 


2.2 


eV 



TABLE II: Constraints on neutrino mass for different exten- 
sions to parameter space for the fiducial CMB+LSS dataset. 



3. AD + all isocurvature modes 

To generate the neutrino velocity mode, a spatially 
varying relative velocity is needed between the neutrinos 
and photons but must be constructed initially such that 
the overall perturbation in total momentum density is 
41| . Such modes therefore require even more exotic 



zero 

generation mechanisms and are met with more skepti- 
cism. Adding in the neutrino velocity mode and permit- 
ting all 4 distinct isocurvature modes to mix, leads to 
a slightly lower neutrino mass constraint of 2.2 eV. The 
distribution of the galaxy bias as shown in figure ^ ex- 
hibits a shift to higher values than in the pure adiabatic 
case. Since the matter power spectrum is dominated by 
the adiabatic mode [27| . including an isocurvature mode 
reduces the overall power P{k), requiring a higher biasing 
of the galaxy spectrum to that of the underlying matter 
distribution. But because our value for b agrees with the 
1.3 ±0.2 obtained for all 4 correlated isocurvature modes 
in the presence of massless neutrinos [45| . we infer that 
the bias is not sufficiently high to require suppression of 
the matter power by such high mass neutrinos. A lower 
value of than that reported in [4oT | must then compen- 
sate for the higher allowed value of neutrino mass. This 
is a strong indication of a degeneracy between ^2 m„ and 
a combination of isocurvature parameters. Note that the 
neutrino mass distribution in figure {2} seems to exclude 
the model with ^m„ = 0. However this feature can be 
attributed to a larger prior space than the massless case 
in which NIV and NID modes are 0. In other words, 
there are effectively more ways of fitting the data with 
massive neutrinos. 




0.0 0.6 1.1 0.0 0.6 1.1 1.6 

I m„ I m„ 

FIG. 3: Bottom panel: The 95 % constraints in the (^m„, 
b) plane for different Gaussian priors on the bias. On left, fib 
is 1.1 (solid), 1.0 (dotted) and 0.92 (dashed). On the right, 
a b is 0.11 (solid), 0.0825 (dotted) and 0.055 (dashed). Top 
panel: The ID marginalized distribution for for the 

same cases. 

IV. THE PATH TO TIGHTER CONSTRAINTS 



If we consider the degradation of the neutrino mass 
constraints which results as the hypothesis space en- 
larges, it seems important to determine how we can im- 
prove our knowledge of ^ m v in the presence of these ad- 
ditional parameters. We have identified two high strong 
degeneracies between m v an d and b which be- 
come important when seeking stronger constraints. Any 
supplementary information that reduces the freedom in 
these degenerate parameters will impact on the neutrino 
constraints. This information must come from indepen- 
dent experiments which measure phenomena on which 
the neutrino mass have negligible effect. 

Stringent neutrino mass constraints can be arrived at 
by limiting the allowable values of b. The upper bound of 
J] m v < 3 eV (95 % CL found by Hannestad [29[ using 
the WMAP lst-year [H and 2dFGRS datasets leav- 
ing b and the normalization as free parameters was low- 
ered to 0.7 eV by the WMAP team using the same data 
by assuming b ~ 1, based on the bi-spectrum analysis of 
the 2dFGRS [H[. Using most recent CMB+LSS data, 
we update this bound to 0.63 eV. The effects of intro- 
ducing different measurements of the bias are illustrated 
in table (jllip . We wish to dissect the effect of including 
bias measurements on neutrino mass constraints in order 
to establish how they need to be improved to better pin 
down ^2 m„ . 

When the analysis is performed using the SDSS data 
its natural distribution is found to be normal with mean 
b = 1.09 and width <jb = 0.11. If we impose different 
Gaussian priors on b keeping the mean fixed and reducing 
the standard deviation to 0.75<7fc and 0.5crfc , the result- 
ing ID ^2 fnu distributions tighten slightly. Because the 
peaks shift towards higher values as shown in figure ([3]) 
the upper bounds at 95% CL are not lowered one might 
expect. With reference to allowed regions in {^2m v , b) 
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Dataset 


E mu for 7V„ = 3.04 
95% Confidence Level (CL) 


CMB 


2.0 eV 


CMB +LSS 


1.3 eV 


CMB+LSS+2dF bias [25] 


1.05 eV 


CMB+LSS+SDSS bias [5] 


0.79 eV 


CMB+LSS+ b = 1 


0.63 eV 



TABLE III: Constraints on neutrino mass for different in- 
formation about the bias 6 for the standard 8-dimensional 
parameter space. 

space, the effect of the imposed standard deviation of b 
on the mass limit progressively weakens, indicative of a 
threshold value of beyond which it has no effect at 
all. If we now lower the mean values y,b of the priors 
but where o\, = 0.1 x /if,, the allowable region in m v, 
b) space is reduced in a relatively linear fashion. Given 
the correlation between b and E m v > a decrease in /ib 
means that a lower b is also favoured and is equivalent to 
weighting the neutrino mass around a correspondingly 
lower value. This shows that in order to improve our 
knowledge of ^ m v a more precise measurement of the 
mean b from future analyses such as weak lensing, will 
be helpful. 

In light of recent results 49] the assumption that power 
spectrum of the galaxy distribution is linearly biased rel- 
ative to the matter power spectrum may be short-sighted, 
with implications for current neutrino constraints. The 
lower matter density favoured by the SDSS DR5 galaxy 
data could be explained in terms of a scale-dependent 
bias. This has the potential to degrade neutrino con- 
straints because the effect of the bias will no longer be 
confined to the amplitude but will now be felt by the 
shape of the linear power spectrum, making it harder 
to break the degeneracy with m v . In work [l9[ in 
which the authors assess the effect of different datasets 
on the neutrino mass limits, discrepancies were pointed 
out between upper bounds on ^ m v obtained using the 
CMB data and either the LSS measurements from 2DF 
Q or SDSS @ . Normalizing the total matter power spec- 
trum without having to make such assumptions can be 
achieved from weak lensing experiments. The limit on 
the neutrino mass found using the measurement of the 
cluster mass function in [50J and the CMB data was con- 
siderably weaker. 

V. CONCLUSIONS 

In this paper, the dependance of constrains on the neu- 
trino mass using cosmological observations (CMB+LSS) 



on the underlying model has been assessed. Before pro- 
ceeding with the analysis we systematically reviewed 
most decisions that go into such an investigation and are 
satisfied that analysis with current datasets is unaffected 
when different neutrino parameterizations are employed. 
The BAO [13J and SN [llj results are shown to impact 
the mass limits signficantly and are thus incorporated 
with caution. 



When parameter space is extended beyond our cur- 
rent standard cosmological model (ACDM), we rule out 
any degeneracy between neutrino mass and spatial cur- 
vature. We also find that the relationship between dark 
energy equation of state w and ^ m„ depends on the 
datasets used. A decrease in w away from —1 can be 
accommodated by an increase in the neutrino mass or 
an increase in n m which rapidly becomes incompatible 
with the SN data. The inclusion of a running spectral in- 
dex no longer tightens the upper bound on m v if the 
SN measurements are removed from the dataset. Both 
of these mentioned results disagree with former findings 
in analyses which include the SN data and highlight is- 
sues regarding the dependence of the results not only 
on the cosmological model but the cosmological observa- 
tions used (see |l9|). Allowing for tensor perturbations 
is found to have minimal effect, however, the analysis ex- 
cluded the possibility of a running spectral index. The 
substantial degradation of the upper mass bound to 2.2 
eV with more general initial conditions points to a degen- 
eracy with the isocurvature parameters, namely the neu- 
trino density mode. These high-^ m v models are how- 
ever ruled out by the SN constraints on Q m which 
bring the limit down to 1.1 7eV. Constraining the value 
of bias impacts significantly on the upper bound, with 
the most stringent constraint of m v < 0.63 eV being 
obtained when taking b = 1. Assuming a known distri- 
bution for the bias amounts to effectively incorporating a 
data set and all accompanying systematic uncertainties. 
For this reason the same discernment should be exer- 
sized. We find that improvements in the accuracy of the 
measurements of the galaxy bias beyond a point will not 
dramatically aid our constraints on ^ m„. In conclusion 
it has been shown that the cosmological constraint on 
neutrino mass is sensitive to many factors and it is only 
once all assumptions have been evaluated can we regard 
the resulting limit as robust. 
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